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We present an experimental procedure to determine the usefulness of a measurement scheme for 
quantum error correction (QEC). A QEC scheme typically requires the ability to prepare entangled 
states, to carry out multi-qubit measurements, and to perform certain recovery operations conditioned 
on measurement outcomes. As a consequence, the experimental benchmark of a QEC scheme is a 
tall order because it requires the conjuncture of many elementary components. Our scheme opens 
the path to experimental benchmarks of individual components of QEC. Our numerical simulations 
show that certain parity measurements realized in circuit quantum electrodynamics are on the verge 
of being useful for QEC. 



I. INTRODUCTION 

The physical realization of a quantum information pro- 
cessor represents a formidable experimental and theo- 
retical challenge. The inevitable presence of noise and 
imperfections in a physical device causes it to deviate from 
its intended ideal evolution. For any complex information 
processing task comprising a large number of elementary 
operations, these deviations can add up to the point of 
completely invalidating the information content of the 
processor. Experimental efforts are constantly improving 
the accuracy in the control of quantum devices and re- 
ducing the level of noise, but this progress can only make 
it so far, and other methods are required to overcome the 
problem of noise and imperfections. 

Quantum error correction (QEC) enables, in princi- 
ple, the coherent manipulation of quantum information 
for arbitrarily long times, despite the presence of noise 
and imperfections in the physical device p~r[3]. This is 
achieved by encoding the information redundantly in an 
error-correction code. The information in the code is 
delocalized, encoded in the many-body correlations of 
the qubits. As a consequence, it becomes less vulnera- 
ble to local noise. To stabilize the information for long 
times, measurements need to be carried on the redun- 
dant systems to detect the presence of possible errors 
and evacuate entropy from the system. Following these 
measurements, a recovery operation can be executed to 
restore the content of the processor to its error-free state. 

Of course, the building blocks of any error-correcting 
scheme are themselves subject to faults, and can some- 
times introduce more errors in the system than they 
remove. The theory of fault-tolerance quantum computa- 
tion [3H5] demonstrates that, provided that the amount 
of noise in each elementary component is below a certain 
threshold value, it is possible to use QEC to decrease the 
level of residual noise after correction below any desired 
value. 

The exact value of this fault-tolerant error threshold 
is not know, but it can be estimated numerically for 
various schemes (e.g. [9] [10]) and theoretical lower bounds 
can be derived for it (e.g. [HI H2]). However, these 
estimates often require simplifying assumptions and fail 



to include all the physical features and limitations of a 
given quantum computer architecture. Theoretical efforts 
are being deployed to remove these assumptions as much 
as possible; correlated errors [T5], slow measurements 
[13], and nearest-neighbor interactions [HI HH HH] are 
a few examples of limitations that have recently been 
incorporated in theoretical estimates of the threshold. 
It remains very likely that the true threshold deviates 
significantly from these estimates, and thus experimental 
investigations of QEC are crucial. 

A handful of experiments [16 - 21] have demonstrated 
the principles of QEC. Part of the difficulty with the 
implementation of QEC is that it requires many building 
blocks, and the failure of one of these could invalidate the 
entire process. Depending on the details of the scheme, 
QEC typically requires the preparation of ancilary qubits 
in a known state, projective measurements, unitary con- 
trol, classical data processing, etc, and moreover each of 
these components often needs to be realized in a very 
short time and with high fidelity. 

The goal of this Article is to demonstrate that some of 
these QEC building blocks can be tested experimentally 
individually, without necessarily realizing a full QEC cycle. 
In particular, we propose a scheme to experimentally 
decide whether a measurement procedure is sufficiently 
good to perform QEC. We hope that this important 
simplification will stimulate experimental research in that 
area, something that could bring new insights to further 
refine theoretical models. We will pay particular attention 
to the measurement process which, in many solid-state 
qubit architecture, can be slow, weak, and noisy. We 
will in particular consider weak parity measurements and 
show that in some conditions such measurements can be 
successfully used in a QEC scheme to reduce the overall 
amount of noise. For concreteness, we will illustrate our 
methods using superconducting qubits in the circuit QED 
architecture [33J HI], but the ideas extend to other devices. 

Another problem that we address and that is often dis- 
regarded is the need for ultra-fast classical data processing 
to assist the QEC scheme. For instance, the continuous 
error-correction schemes of Ref. [24 is very effective at 
suppressing errors in theory, but the data processing it 
entails would require a speed-up of current classical data 
processing of several orders of magnitudes before becom- 
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ing practical. Fault-tolerant schemes tailored for slow data 
processing [13] . more cleaver data- processing algorithms 
[2"5] . or even fast FPGA electronics could surmount these 
roadblocks, but they nevertheless remain an obstacle to 
current experimental investigations of QEC. We will show 
that for the purpose of benchmarking, all data processing 
can be delayed until after the experiment at no cost. 

This paper is organized as follows. We first present in 
Sec. [IT] an overview of QEC and highlight some difficul- 
ties associated to implementation of these ideas with the 
measurements typically realized in solid-state qubits. In 
Sec. we propose a figure of merit that isolates the use- 
fulness of a measurement scheme in a QEC protocol, and 
propose an experimental procedure to evaluate it. The 
validity of this experimental procedure is demonstrated 
in the appendix. Section |IV| presents an numerical ap- 
plication of our method to a recently proposed parity 
measurement in circuit quantum electrodynamics archi- 
tecture [26]. 



II. QUANTUM ERROR CORRECTION 
A. General description 

A QEC cycle usually breaks into three steps, the last 
one being to some extent optional. In a first step, the 
system is initialized to a code state S C. The quantum 
code C is a subspace of the Hilbert space of n physical 
qubits that can be specified by a set of operators S a , called 
stabilizers. The stabilizers represent constraints on the 
encoded information in the sense that is a code state if 
and only if it gives the outcome +1 when any stabilizer S a 
is measured, i.e. C = : S a \i/j) = \ip)}- The dimension 
of C determines the number of encoded qubits. When 
dim(C) — 2 k , we say that the code encodes k qubits, and 
that it has a rate k/n, the ratio of encoded qubits per 
physical qubits. For simplicity, we will here limit our 
discussion to codes encoding a single qubit, k = 1. 

After the information is encoded, the physical qubits 
are subject to a noisy evolution. The second step of QEC 
consists of verifying that the state remains inside the code 
after this evolution. This is achieved by measuring the 
stabilizers S a , a measurement outcome that differs from 
+1 signalling the presence of errors. The operators S a are 
usually elements of the Pauli group — i.e., constructed 
from tensor products of the three Pauli matrices o~ x , o~ y , 
and a z and the identity — and can therefore only yield 
±1 measurement outcomes. The collection of these ±1 
measurement outcomes, one per stabilizer S a , is called 
the error syndrome. 

The third and last step consists of reverting, to the 
greatest possible extent, the effect of the noise. Such a 
recovery operation can be inferred from a prior statisti- 
cal description of the noise model — as given, e.g., by a 
relaxation time Ti and a dephasing time T v — and the 
additional information obtained from the error syndrome. 
It is common to assume that the errors E are also ele- 



ments of the Pauli group, so the statistical evolution of 
the system is captured by the quantum map £ 

E 

where Pr(E) represents the probability that the Pauli 
error E = E> has affected the system. In that case, the 
recovery operation could consist of applying to the system 
the Pauli operator E that has the largest probability 
Pr(-E) and which is consistent with the observed error 
syndrome. This will be effective since all elements of the 
Pauli group square to the identity, so applying the same 
error twice has no net effect on the system. In certain 
circumstances, this third step can be partly omitted; it 
is often sufficient to keep in mind that the information 
store in the encoded qubits has been modified by some 
known unitary transformation E — i.e. perform a basis 
change — and unnecessary to revert its action. This is 
sometimes referred to has changing Pauli frame [TTH |2"T] . 



B. Difficulties in application to solid-state qubits 

In practice, many of the usual assumptions described 
above are not met in solid-state devices, and this leads 
to additional burdens for experimental investigations of 
QEC. First and foremost, breaking the QEC cycle into 
discrete steps — a noisy evolution followed by a projective 
measurement — is not always possible. Measurements of 
solid-state qubits are often slow, and errors will occur on 
the system during the course of these measurements. 

In addition to being slow, measurements are noisy. For 
example, often, the experimental outcome for the mea- 
surement of a Pauli operator will be a continuous and 
noisy signal rather than a binary ±1 value. This is, for 
example, the case with an heterodyne measurement in 
circuit QED [28]. Of course, the noisy signal can be inte- 
grated over a period of time and discretized to be treated 
as binary, but this would entail a loss of precious informa- 
tion about the errors, as we will see below. Measurements 
can be omitted, but at the price of performing multi-qubit 
gates (outside the Clifford group) and resetting auxiliary 
qubits [2"T] . 

Lastly, the noise affecting solid-state devices (and other 
qubit architectures) does not correspond to a simple statis- 
tical distribution Pr(E) over elements of the Pauli group. 
Instead, we usually model the evolution of these systems 
by a master equation of the Lindblad form. This more 
complicated evolution has several consequences. In gen- 
eral, it requires an error-diagnostic protocol far more 
complex than the one described above for Pauli operators 
that simply consisted of optimizing Pr(_E) over a subset 
of Pauli operators. In addition, the optimal recovery op- 
eration found by the error-diagnostic protocol will not 
be a simple unitary transformation, but instead will be 
a general quantum map. The physical realization of this 
map could be beyond today's technologies. Moreover, 
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because the recovery will typically not be a unitary trans- 
formation, it cannot be regarded as a basis change and 
it becomes imperative to implement it physically, an ad- 
ditional experimental burden. One can always pretend 
that the noise model is a Pauli channel and perform error 
correction following the usual approach. However, this 
will inevitably entail a loss of performances, so it does not 
well reflect the usefulness of the measurement scheme. 



C. Two examples of quantum codes 

We will illustrate our methods to test QEC building 
blocks using two different error correcting codes. The first 
one is the bit flip code [T] . It uses n = 3 qubits and has 
stabilizers Si — a z ® o z ® I and S2 = I ® cr z ® a z . Thus, 
the code C is spanned by the two states |0) = |000) and 
|1) = 1 1 1 1) . As its name suggests, this code is effective at 
correcting bit flip errors that invert the two basis states 
|0) <-> corresponding to the quantum operation a x . 
On an encoded state a|0) + /3|1), flipping the first qubit 
would lead to a state of the form a|100) + /3|011), with 
error syndrome S\ = — 1 and 5 2 = +1. Similarly, flipping 
the second qubit produces the syndrome (—1,-1) and 
flipping the third qubit produces (+1,-1). Thus, any 
single-qubit flip can be uniquely identified and corrected 
by an additional flip. 

Flipping two qubits leads to the same syndrome as 
flipping only the complementary qubit. For instance, an 
error that flips the first two qubits E = a x <g> a x <g> I 
leads to the syndrome (+1, —1), and would be mistaken 
for the error E' = I ® I ® <j x flipping only the third 
qubit. Correcting the error E with the misdiagnosed 
flip E' would have the net effect of flipping the encoded 
state from a\0) + to a\l) + /3|0). Thus, the code can 
correct any one-bit flip, but fails at correcting two-bit 
flip. For a noise model where each qubit is independently 
flipped with probability p, the code would have a failure 

probability i^)p 2 , a net gain when p < ( 2 ) = |. 

In practice however, bit-flip is not the dominant source 
of noise in solid-state devices, so the bit-flip code is not 
very relevant. We nonetheless choose to use this code to 
illustrate some of our schemes because it is conceptually 
very simple. Note that in this paper, whenever we use 
this code, we will artificially change the noise model of 
the device, adapting it for the purpose of illustration with 
the bit flip code. 

The second code we study is tailored to deal with 
relaxation processes [29j , which are more directly relevant 
for solid-state qubits. Indeed, relaxation is the main error 
channel in superconducting qubits, and pure dephasing 
noise is often negligible in comparison |30j . The code uses 
n = 4 qubits and has stabilizers S% — a z ® o z ® I ® I ', 
S 2 = I ® I <S> <J Z ® a z , and S 3 — a x ® a x ® a x ® a x . 
Thus, the code C is spanned by the two basis states 
|0) = ^(|0000) + 11111)) and |1) = ^(|0011) + |1100». 
Although the workings of this code is more intricate, one 
can gain some intuition by noting that when the first qubit 



relaxes, corresponding to the operator a ® I ® I <g> I, the 
initial state a|0) + /3|1) maps to 

a|0111) +/3|0100), (2) 

retaining all the information about the encoded state. 



III. BENCHMARKING THE MEASUREMENT 
A. Isolating the measurement 

Stabilizers can always be measured indirectly by per- 
forming an appropriate multi-qubit gate, followed by 
single qubit measurement. Alternatively, they can be 
measured directly using some sort of generalized parity 
measurement. For instance, the parity measurement pro- 
posed in [2H l3"lTf3"3"] can be used to measure two-qubit 
stabilizers, which includes all the stabilizers from the 
examples described above with the exception of S3 in 
the last example. While both approaches are viable, our 
scheme will focus on such direct stabilizer measurements. 

The question we would like to answer now is: Can a 
given physical realization of a parity measurements be 
successfully used in a QEC scheme to reduce the overall 
amount of noise? One way to answer this question is 
to use the parity measurement in an experimental QEC 
protocol and check that the resulting noise rate is re- 
duced. However, this method would provide an answer 
that depends on all the components of the QEC proto- 
col. For instance, perfect parity measurements combined 
with noisy recovery operations could lead to an overall 
useless scheme. Here we will show how to decouple these 
questions and propose a method to evaluate the quality 
of the measurement alone. 

To decouple the measurement process from the other 
components of the QEC scheme, it is useful to consider 
the following two gedanken experiments. In the first ex- 
periment, the system is prepared in a random unknown 
code state \tp), and subject to noisy evolution for a time 
T, resulting is a mixed state p — £ ('I'). Here and through- 
out, we use capital Greek letters to denote the density 
matrix corresponding to a pure state \t = £ is the 

completely positive trace-preserving (CPTP) map that de- 
scribes the evolution of the system for time T. After this 
time, the optimal recovery operation is performed on p to 
maximize the fidelity to the initial state. This recovery 
operation can be anything allowed by the laws of quantum 
mechanics, including measurements, unitary transforma- 
tions, coupling to an external system, etc. Thus, it can 
be described by a CPTP map 1Z. In other words, one 
chooses the CPTP map 1Z that maximizes the quantity 
F R (ip) = (ip\K(p)\ilj) = (il;\7i(£{^))\^) : averaged over all 
code states \ip). The average performance of this scheme 
can be summarized by the average fidelity 

Fx =max [ Fu^dijj. (3) 
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In the second thought experiment, the system is again 
prepared in a random unknown code state \ip). Then, con- 
tinuous weak parity measurements of the stabilizers {Si} 
arc performed for a time T to stabilize the encoded states. 
The outcome of these measurements is a collection of 
continuous current J(t) = { Ji(t)}, one for each measured 
stabilizer generator. Of course, while these measurements 
are being performed, the system undergoes its noisy evo- 
lution and suffers from measurement back-action. The 
combined effect of the noisy evolution and the measure- 
ments results in a state p = £j(^), where £j is a com- 
pletely positive (CP) map that depends on the observed 
currents J . At time T, after the evolution, one can apply 
a recovery 1Z which yields a fidelity to the initial state 
Fn(il>,J) = (mj(p)W = WKj(£j(*)M. Clearly, 
the optimal choice of 1Z will depend on the observed cur- 
rents J . The average performance of this scheme can be 
summarized by the average fidelity 



F 9 =E 



max 
n 



F n (^,J)P(iP\J)dip 



(4) 



J 



where E[-]j stands for the average over the output cur- 
rents J, whose distribution depends on the initial state 
) and can be modeled by the conditional probabil- 
ity P{J\iji). The probability P(ip\J) is obtained from 
P{J\ijj) via Bayes rule. We denote TZj the recovery that 
maximizes the average fidelity given an output current 
J. 

Our original question of whether the parity measure- 
ment can be used in a QEC scheme to reduce the over- 
all amount of noise can then be answered positively if 
F\ < F2- Indeed, the only difference between Eq. ([3| and 
Eq. Q is that in the latter, the system's dynamics is 
altered by the presence of measurements and one uses 
the outcome of these measurements to design the optimal 
recovery. The presence of the measurement can lower the 
fidelity because of the measurement back-action, but it 
can also increase the fidelity because it brings informa- 
tion about the stochastic evolution of the system. Thus, 
F\ < F% means that including the measurement had an 
overall positive effect. 



B. Optimal recovery 



state |0) = ^(|0) ® |0)r + |1) O \\) R ) that is maximally 
entangled between the code space and some auxiliary ref- 
erence qubit. Here, the subscript R denotes the reference 
qubit. The entanglement fidelity describes how well this 
entangled state was preserved by the combined effect of 
the noise £j and the recovery TZj: 



F e {J) = {cj>\Kj{£j{m<i>) 



(5) 



where it is understood that both 1Z j and £ j act only on 
the system, i.e. the reference qubit evolves trivially. The 
average fidelity and the entanglement fidelity are related 
by F = where d is the dimension of the code space, 

in the present case d ~ 2. 

Thus, maximizing the average fidelity is equivalent to 
maximizing the entanglement fidelity. This is an impor- 
tant conceptual simplification because it only requires 
optimizing the fidelity for one given initial state — an 
entangled state between the system and the reference 
qubit — and therefore leads to a simple numerical proce- 
dure. All we need to do is to numerically integrate the 
master question for the system including noise, given the 
initial state \<j)) and taking into account the observed cur- 
rents J [26l [34] . The outcome of such a simulation is 
a mixed state which we will denote = (£j ® /)($). 
For a recovery operation 7Z, the entanglement fidelity can 
now be expressed as 



F e = (cf>\(K®I)(n £ j)\cl>) 



(6) 



where Q,n = 7Z($>). Thus, the optimal recovery can be 
found by maximizing Tr(fi K t^£j) over fi^, subject to 
the constraint that 1Z is a CPTP map, i.e. FIr > 
and TtrQh — I, where Tr# stands for the partial trace 
over the reference qubit. This optimization problem is a 
semi-definite program [35], and can be efficiently solved. 

To summarize, the optimal recovery TZj for a given set 
of output currents J7, can be found by first integrating the 
evolution equation for an input state entangled between 
the code space and a reference qubit, and then solving 
a semi-definite program. Both of these can be realized 
in polynomial time as a function of the dimension of the 
Hilbert space of the system and the duration T of the 
experiment. 



In this section, we explain how to mathematically con- 
struct the optimal recovery maps 1Z needed for the two 
gedanken experiments described above. We will focus 
on the second type of experiment, corresponding to F2 
and where continuous parity measurements are performed. 
The case without measurements follows trivially by re- 
moving all references to the output currents J . The 
average fidelity to the initial state F is a good figure of 
merit because it has a simple interpretation. However, 
it is mathematically more convenient to use a slightly 
different figure of merit called the entanglement fidelity 
F e defined as follows. Imagine preparing the system in a 



C. Delayed tomography 

So far, we have described how to compute the optimal 
recovery operation 1Z to revert the effect of some evolution 
£ , that may or may not contain contributions from parity 
measurements. This recovery is chosen to maximize the 
average fidelity over the code space. The usual way to 
assess this protocol experimentally is to perform quantum 
process tomography on the map ¥\Rj£j]j (or 1Z£ in the 
absence of measurements), describing the joint effect of 
the measurement and the recovery. This can be realized 
experimentally as follows. 
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1. Prepare a code state \4>i)- 

2. Let the system evolve according to £j. 

3. Compute the optimal recovery IZj. 

4. Apply the optimal recovery TZj to the system. 

5. Measure some observable Aj. 

Repeating this sequence many times provides an estimate 
of E[Tr{A J -7ej-(£j-(|^XV ) i|))}]j'- Varying over the initial 
state \ipi) and final measurement Aj provides complete 
information about the map K[1Zj£j]j. From this com- 
plete description, on can directly compute the average 
fidelity. 

There are two obvious problems with this approach. 
First, the optimal recovery IZj might be beyond what 
can be realized with today's technologies, so step 4 cannot 
be realized. Second, when £j is the evolution resulting 
from the combined effects of noise and parity measure- 
ments, the optimal recovery IZj depends on the observed 
measurement currents J. In that case, the procedure to 



This procedure is repeated M times, and the average 
entanglement fidelity F e — ¥\F e {J)]j can be estimated 
by 



compute TZ j described in Sec. Ill B must be realized in 
real time. With the integration of the evolution equation 
of a 3-4 qubit system taking several seconds, the system 
will have long completely decohered before the recovery 
operation can even be applied. 

We propose a method to directly estimate the average 
fidelity that circumvents these problems. The intuitive 
idea is to invert steps 4 and 5 in the above procedure, i.e. 
measure a complete set of observables before applying the 
recovery. So crudely, we can imagine that we are char- 
acterizing the state £j(&) before applying the recovery, 
and then computing what would have been the resulting 
fidelity had we applied the recovery. The problem with 
this crude explanation is that the state £ j($) changes at 
every run of the experiment because it depends on the 
observed currents Jj(t). Thus, this naive approach only 
reveals the average state E[£j(<i>)] i 7, which washes away 
the information content of Jj(t) that should have been 
used to determine the optimal recovery. 

Nevertheless, the essence of the idea is right and, bor- 
rowing some ideas from Refs. [36l[37j, we will show that 
the following scheme can be used to estimate the average 
fidelity. 

1. Choose a randomly among the Pauli operators and 
a random r = ±. Prepare the n physical qubits in 
the state |V>J), that is the eigenstate with eigenvalue 
r of the encoded operators <r, i.e. a state among the 
set{|0),|l),|0)±|T),|0)± l |l)}. 

2. Let the system evolve according to £ j, and collect 
the currents J . 

3. Measure the observable chosen uniformly at 
random from the set of n-qubit Pauli operators, 
obtaining the result v. 

4. Compute the optimal recovery operation IZj. 
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2^AMTV(n; t Bfl®aW), (7) 



where the superscript [£] makes reference to the £th re- 
alization of the procedure. In this altered scheme, the 
experimental work ends at step 3; step 4 and the compu- 
tation of fidelity can be delayed as post processing. The 
validity of Eq. Q is demonstrated in appendix [A] 



IV. EXAMPLE: PARITY MEASUREMENT IN 
CIRCUIT QED 

A. Parity measurement 

As a concrete example of QEC measurements bench- 
mark, we focus on the 3 and 4-qubit codes discussed above 
and consider a realization in circuit QED. These ideas 
are however general and apply to other systems. We will 
consider the continuous two-qubit parity measurement 
proposed in Ref. for circuit QED. Note that two-qubit 
parity measurements are sufficient to implement the bit- 
flip code, but insufficient for the relaxation code that also 
contains the stabilizer S3 = cr x ®o- x ®a x ®a x . In our case 
study, this third stabilizer will simply be ignored, only the 
stabilizers that involve a simple two-qubit parity being 
measured. Since S3 is responsible for the correction of 
phase errors that occur on a longer timescale T 2 — 2Xi , we 
expect that measurement of the two-qubit parities alone 
can lead to an improved error rate, and this is confirmed 
by our numerical simulation. 

Before presenting the numerical results, we first give 
a brief overview of circuit QED. In this architecture, a 
superconducting qubit is dipole coupled, with strength 
g, to the zero-point electric field of a transmission-line 
resonator [2"2l |2"3"] . In the regime where the detuning A 
between the qubit transition frequency u a and the res- 
onator frequency uj r is large with respect to the coupling 
strenght, |A| = \ui a — u r \ >• 3, the electric-dipole interac- 
tion leads to a dispersive coupling of the qubit and the 
resonator of the form x a ^ ao ~z, where x = 9 2 /^. This 
interaction can be interpreted as a qubit-state depend 
shift of the resonator frequency to uj r + x°~z- As a result, 
a coherent tone of frequency uj m ~ ui r and amplitude 
e m driving the resonator will displace the intra-resonator 
field, initially assumed to be in the vacuum state, to a 
qubit-state-dependent coherent state \a a ) of amplitude 
a a = —e m / [u! r + (— 1) CT X — w m ~ in/2]. In this expres- 
sion, k is the resonator photon-loss rate and a = {0, 1} 
stands for the two qubit states. These two coherent state 
can be resolved by homodyne detection of the coherent 
drive transmitted or reflected by the resonator, thereby 
realizing a measurement of the qubit [38] . 
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FIG. 1: (color online) Phase space illustration of the stationary 
states \a aa i) for: g\ = —g% = —15k and Xj ~ 1.5k. The 
measurement drive is such that cj m = lo t and e m = re/2. 



A similar situation holds if multiple qubits are dis- 
persively interacting with the resonator. For exam- 
ple, in the presence of two qubits, the resonator fre- 
quency is shifted to the two-qubit state dependent value 
L °r+Xi< J z ( 3l+X'2l ( 3< J z, where \j is the dispersive coupling 
of qubit j. This leads to four possible qubit-state depen- 
dent coherent states \a a ), where now a — {00,01, 10, 11}, 
and that can be resolved in a homodyne detection |39j . 

As shown in Refs. (2^1110], it is possible to adjust the 
frequency of the measurement tone, the dispersive cou- 
plings and the frequency and phase of the local oscillator 
used in the homodyne detection such that the measure- 
ment only reveals information about the parity of the 
joint two-qubit states. This is illustrated in Fig. [I] where 
the parameters have been chosen such that the coherent 
states corresponding to even qubit parity {|00), |11)} over- 
lay, while those corresponding to odd parity {|01), |10)} 
almost do 1 . In this situation, homodyne detection of the 
Q quadrature of the resonator field will yield information 
about the qubit parity, and nothing about the individual 
qubit state. In other words, for the settings of Fig. [TJ 
the measurement operator in circuit QED takes the form 
er z (g)cr z , exactly what is required for the stabilizers S± and 
S2 of the two quantum error correcting codes discussed 
above. Moreover, these two stabilizers can be measured 
simultaneously, for example, by using a two-resonator 
setup where the pair of qubits to be jointly measured is 
fabricated in the same resonator. For the 3-qubit code, 
this means that one of the qubit must be coupled to 
two resonators, something that has already been realized 
experimentally |41j . 



1 It is not possible for the coherent states corresponding to odd 
parity to perfectly overlay if the even ones do, and vice-versa. 
With the parameters chosen for Fig.JT] this results in measurement- 
induced dephasing of superpositions of the odd parity states |26| . 
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FIG. 2: (color online) a-b) Fidelity vs time for optimal QEC 
with (blue line) and without (green line) parity measurements, 
textbook QEC (purple line) and unencoded single qubit evo- 
lution (yellow line). Panel a) is a zoom of panel b) in the 
region close to unit fidelity, c) Fidelity vs detection efficiency 
r\ taken at t = 48 ns as indicated by the vertical dashed line in 
panels a) and b). In all cases, the parameters are x/^ n = 120 
MHz, k/27t = 50 MHz, ~/ x /2tt = 5 MHz. When measurement 
is present, 77 = 1 [in panels a) and b)] and e/2-K = 40 MHz. 
The results where averaged over 8000 trajectories. 



B. Numerical results 

Fig. [2] presents the fidelity as a function of time for the 
three-qubit code under continuous measurement of the 
syndromes Si and S%. These results were obtained by 
numerical integration of the stochastic master equation 
(SME) found in Eq. (2) of Ref. and using the param- 
eters indicated in the caption. This SME describes the 
evolution of two qubits coupled to the same resonator 
and under continuous parity measurement realized, as 
discussed above, by homodyne detection of the signal 
transmitted through the resonator. Since the 3-qubit 
code protects against bit flips, and not relaxation, we 
replaced the qubit damping term in this SME by a term 
representing symmetric bit flips. We have also dropped 
pure dephasing since both the 3 and 4 qubit codes do 
not protect against this noise source (and many super- 
conducting qubits in circuit QED suffer only very weakly 
from pure phase damping [SU])- These are the only modi- 
fications made to the results of Ref. 26J . 

In this Figure, the blue and green lines correspond 
to the fidelity after optimal recovery, respectively in the 
presence and absence of continuous weak parity measure- 
ments. These results are compared to the fidelity of a 
single unencoded qubit (yellow line). For the chosen pa- 
rameters, we see that it is always beneficial to encode 
the information (blue and green line are above the yel- 
low one), but that the benefit of weak measurements 
depends on the duration of the experiment. Thus, even 
with optimal recovery, continuous syndrome measurement 
does not always help. Indeed, below ~ 25 ns, measure- 
ment backaction [35] disrupts the encoded state more that 
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to an unencoded qubit, the green line to optimal recovery 
in the absence of continuous measurements and, finally, 
the blue line to optimal recovery with continuous measure- 
ment of Si and 52- Compared to the three-qubit code, 
the case with measurement always do as good, or better, 
than without measurement. This is because measurement- 
induced dephasing acts in the odd-parity subspace |26j . 
This, combined with the fact that the four-qubit code- 
words are of even parity and that relaxation never bring 
a superposition of encoded states to a superposition of 
odd-parity state 2 , yields the enhanced result compare to 
the three-qubit code. 



C. Possible experimental realization 



FIG. 3: (color online) a) Fidelity vs time for optimal QEC 
with (blue line) and without (green line) parity measurements, 
and unencoded single qubit case (yellow line), b) Fidelity vs 
detection efficiency rj taken at t = 48 ns as indicated by the 
vertical dashed line in panel a). The parameters are the same 
as in Fig. [2] 



the acquired information helps. Beyond that time, the 
measurement record can be used to improve the optimal 
recovery. 

The purple line shows the fidelity obtained by a direct 
application of textbook QEC. The continuous weak stabi- 
lizer measurements are performed for a time t, after which 
the output signal are integrated and round-up to a ±1 
value. The simplest Pauli operator consistent with this 
round-up error syndrome is applied to the state. At early 
times, the fidelity of this scheme approaches 1 /4. This is 
because the integration of the output signal is essentially 
random, so the recovery will essentially consist in flipping 
one of the 3 qubits or none at all, each option having 
probability 1/4. Since at early times it is likely that the 
qubit has not suffered any errors at all, only 1 out of 
these 4 options is correct, explaining the 1/4 fidelity. The 
fidelity reaches a maximum at a subsequent time, when 
the syndrome measurements reveal useful information 
about the error. Yet, it does not reach the fidelity under 
optimal recovery. This is because, as explained above, 
the optimal recovery can be any CPTP map, while this 
textbook application is limited to a recovery that flips 
one of the three qubits. As intuitively expected, the time 
scale where this maximum occurs corresponds roughly 
to the time scale where weak measurements provide an 
advantage under optimal recovery; prior to this time the 
measurement signal is essentially random. 

Fig. [3] shows similar results for the four qubit code. 
Since this code protects against relaxation (and not sym- 
metric bit flips), this is the type of noise that was numer- 
ically simulated to generate these results. This situation 
is therefore closer to the experimental reality. Using the 
same color code as in Fig. [2j the yellow line corresponds 



Is this protocol experimentally realistic? First, we 
emphasize that, with delayed tomography, the recovery 
operation need not be applied in the laboratory. It is 
only applied, after the fact, on the tomographically re- 
constructed state. The most stringent requirement is on 
the quality of the parity measurements, something that 
depends on the details of the implementation. In this 
respect, all the circuit QED parameters that we have 
chosen are within experimental reach, except one: the 
detection efficiency r\. This parameter corresponds to 
the the efficiency with which photons used for the syn- 
drome measurements are detected. The results presented 
in Figs. [2] and [3] have been obtained by setting 77 = 1, 
above what can now be realized. Indeed, most exper- 
iments use cold amplifiers with a rather low detection 
efficiency r\ f=a 0.05. Recent experiments with near quan- 
tum limited amplifiers have reached rj « 1/3 [42] . and 
values approaching unity should be realizable in the near 
future. To evaluate the importance of this parameter, 
the fidelity, evaluated at t = 48 ns, is plotted versus 77 in 
Figs. [2]:) and (3)3). For the three qubit code, because of 
the measurement-induced dephasing, the weak syndrome 
measurements do not help until a rather high value of 
rj ~ 0.85. For the more experimentally relevant case of 
the four-qubit code, because this backaction does not 
damage the encoded state, measurement helps even at 
low efficiency. 

The two other parameters worth commenting on are the 
coupling strength g and the relaxation rate 71 . Our choice 
of x/2tt = 120 MHz and A = g/A = 1/10 3 correspond to 
a qubit-resonator coupling strength of g/27r = 1200 MHz, 



2 Note for exampfe that the state in Eq. can be written as 
|01}(S'(a|ll}+/3|00}), and contains no superpositions of odd parity 
states that would undergo measurement-induced dephasing. 

3 The value of A = g/A = 0.1 is chosen such that the system 
is in the dispersive regime where the results of Ref. I26| hold. 
For the same reason, we have also chosen the measurement tone 
amplitude e such that the average photon number in the resonator 
is always smaller than about a tenth of the critical photon number 
n crit = (A/2g) 2 [H]. 
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a value that is beyond current experiments by about a 
factor of four [21]. An architecture to realize this type of 
'ultra-strong' coupling has been theoretically proposed [35] 
and tested experimentally with g/2ir ~ 600MHz [331 145] . 
The value used here is thus not beyond reach. Finally, the 
relaxation rate 71 /2ir = 5 MHz used here is more than an 
order of magnitude larger than recently achieved [35J [37J . 
Adjusting this rate to the best current experimental value 
would relax the constraints on all the other parameters. 

Lastly, our scheme is sensitive to state initialization 
errors. These are estimated to be at the percent level [48] 
but can be reduced by using measurements to initialize 
the qubits [351150]. 



V. DISCUSSION 



on the error-corrected qubit (but correlated instead to 
error history on the ancillary qubit). While some of the 
techniques we have proposed could be extended to an 
indirect measurement setting, it is not clear how to decou- 
ple the measurement imperfection from the imperfections 
arising in the unitary control in a measurement scheme 
that relies on unitary control. 4 

Finally, another interesting error correcting code to 
study with this architecture is the Bacon-Shor code [5TT - 
[53] . The main motivations for this code is that it is 
spatially local in a 2D architecture, requires only two-qubit 
parity measurements, and can correct all type of errors. 
The complications arise because the parity measurements 
need to be in complementary basis, and our theoretical 
treatment would need to be generalized to subsystem 
codes [5T] [53]. 



The method we have presented offers the possibility 
to experimentally determine if a measurement scheme 
is useful for a given QEC protocol, without the need to 
realize the full QEC protocol. In particular, there is no 
need in our method to implement the recovery operation 
on the system and all data analysis are delayed until after 
the experiment — in particular there is no operation that 
is realized on the system conditioned on the outcome of 
previous measurements. 

There are obvious limitations to our scheme. On the 
one hand, a negative answer does not imply that the mea- 
surement scheme is of no use in QEC, but merely as used 
in a prescribed scheme. In principle, there can exist QEC 
schemes that are more or less sensitive to measurement 
imperfections than other. On the other hand, a positive 
answer only indicates that the measurement is useful in 
principle, but other limitations, such as imperfect unitary 
control and slow data processing, could render the en- 
tire scheme useless. Despite these limitations, we believe 
that our scheme is of interest at this early stage of infor- 
mation processing devices, and in particular for devices 
where distinct measurement schemes are currently being 
developped. 

The type of syndrome measurement that we have con- 
sidered in this Article are of a special type that we call 
direct: the output signal at time t is directly correlated 
with the value of the measured observable at that time 
(or a slightly earlier time) . This type of measurement is 
to be distinguish from an indirect measurement where 
an ancillary system is coupled to the measured qubits 
and measured at a later time. For instance, a syndrome 
measurement can be realized by applying cnot gates to an 
ancillary qubit and then measuring that ancillary qubit 
in the computational basis. In a direct measurement, 
the signal history correlates with the noise history, so 
there is information to be gained by analyzing the entire 
signal rather than simply considering its average. In an 
indirect measurement, all the correlations between the 
error-corrected qubits and the ancillary qubit are estab- 
lished at the beginning (and on a very short time scale), 
so the signal history is uncorrelated to the noise history 



Acknowledgments 

We thank K. Lalumiere, M. P. da Silva and M. Bois- 
sonneault for discussions and help with the numerical 
calculations. This work was partially funded by FQRNT, 
NSERC, the Alfred P. Sloan Foundation, and CIFAR. 



Appendix A: Estimation of the average 
entanglement fidelity 

For a complete set of observables Aj satisfying 
Tr(AjAk) = ddjk, for example Pauli operators where 
d = 2, we can express the entanglement fidelity as 



F e = Tr(ft K t fi £ 



(Al) 



Let us denote by \ipl) E {|0), |1), |0) ± |1), |0) ± the 
eigenstates of the Pauli operators, i.e. cr\i/j^} = t|V£) and 
similarly, let us denote by |V£) £ {|0), |T), |0) ± |T), |0) ± 
the encoded versions. We can decompose the Pauli 
operator Aj = B k ® a into a component B k on the n 
physical qubits forming the QEC code and a Pauli matrix 
a acting on the reference qubit. This enables us to write 

Trifled) = Tr{B k £j(*+)} - Tr{B k £j(^ -)}. (A2) 



Combining with Eq. (Al), we obtain 



F, 



EE E rTv{n n , B k ®a)Tv{B k £j{^l)}. (A3) 



a r— ± 



4 In the circuit QED architecture we have considered, the cavity is 
coupled to the qubit throughout the measurement process. Even 
though the cavity serves as an ancillary system, the measurement 
cannot be decomposed into a two stage process (unitary coupling 
followed by a measurement) and is thus of the direct type. 
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Similarly, we can take a spectral decomposition of B k ■■ 
P/l — to obtain 

F - = E E E rvTr(il n , B k ® a)Tr{P^ J (*;)} 



cr t,i/=± 



(A4) 

The first term TrffZ—t ® cr) of this expression can 

j 

be evaluated numerically given knowledge oilZj. Let us 
denote 



Pr(^, t, fc, <t) = 



Tr{P fc ^(^)} 
2 • 4"+! 



(A5) 



random among the 4™ n-qubit Pauli operators; choosing 
a uniformly at random among the 4 single-qubit Pauli 
operators; r uniformly at random between ±1; and v 
at random according to the probability Tr{P^fj(4 , J)}. 
These probabilities correspond exactly to probability with 
which the experiment described in Sec. |III C| generates 
the values k, a, t, and v. 

Finally, rewriting 



This can be interpreted as a probability distribution; it 
is positive and sums to 1. Moreover, one can efficiently 
sample from this distribution by choosing k uniformly at 



we recognize Eq. as a Monte Carlo estimate of the 
entanglement fidelity F e . 
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